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We consider a model for population dynamics such as for the evolution of bacterial colonies 
which is of the Fisher type but where the competitive interaction among individuals is non-local, 
and show that spatial structures with interesting features emerge. These features depend on the 
nature of the competitive interaction as well as on its range, specifically on the presence or absence 
of tails in, and the central curvature of, the influence function of the interaction. 



PACS numbers: 87.17.Aa, 87.17.Ee, 87.18. Hf 



Population dynamics covers a wide spectrum of fields. 
The formation of patterns in the evolution of bacterial 
colonies provides an example QI^IHQiI^- Such studies, 
e.g., the growth of an initial nucleus of cells, are impor- 
tant within a clinical framework since the control of bac- 
terial infections might be better performed if the dynam- 
ics of the growth of bacteria were understood at a basic 
level. The Streptococcus pyogenes or Group A strepto- 
coccus, that has been made famous in recent times by 
media reports of flesh-eating bacteria which claim the 
lives of up to 25% of its victims, provides one example. 
Processes related with cancer development provide an- 
other. 

When dealing exclusively with the evolution of cell 
population and its spatio-temporal features, it is usual to 
focus attention on some processes such as reproduction, 
competition for resources, and diffusion, and to neglect 
others such as mutation. The Fisher equation is, there- 
fore, a good starting for such studies. We show below 
new results concerning pattern formation that arise from 
a natural extension of the Fisher equation. This differen- 
tial equation considers, besides the diffusion process with 
coefficient D, growth of the population at rate a and a 
limiting process controlled by 6, associated generally with 
competition or struggle for resources 0|: 
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Our generahzation consists in incorporating non-local 
effects in the competition terms: 
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Here /cr(x, y) is a positive function or distribution which 
we call the influence function, characterized by a range 
a, and normalized in the domain f2 under investigation. 
The physical origin of non-local aspects in the compe- 
tition interaction is easy to understand. For instance. 



in the case of bacteria, the diffusion of nutrients and/or 
the release of toxic substances can cause non-locality in 
the interaction. While non-local competition has been 
mentioned earlier ^ iS] , we show below consequences of 
the non-local terms that have not been reported earlier, 
specifically, striking features of the dependence of the 
patterns on the nature as well as range of the infiuence 
function. 

It is known that no patterns appear in the extreme 
local limit /^(x, y) = 5{x — y) which reduces (01 to CJ. 
No patterns appear in the extreme nonlocal limit either, 
where fcr{x) is a constant, as can be seen from an explicit 
analytical solution given by one of the present authors 0| : 
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where \& (z, t) is the standard Gaussian propagator 
{4:7rDt)~^/^ exp{— / ADt) of the diffusion equation, and 
uo is the integral over all space of the initial density 
uo{x). In the intermediate case, however, patterns do 
appear. The simplest forms for the infiuence function 
are a Gaussian and a square distribution. The latter is 
simply a normalized constant function within a range and 
vanishing outside. The former is given by 

{x-yf 



fa{x,y) = a{x) exp 



(4) 



in the 1-dimensional case. In the case of periodic bound- 
ary conditions with spatial period L, the normalization 
factor a{x) is a constant, ^ = y^cierf whereas, 
for zero-fiux boundary conditions, it is space-dependent: 
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Is is helpful to characterize the influence function by its 
width at the origin. For periodic boundary conditions, in 
which case fa depends on the difference z = x — y, this 

width is given by E = [—d'^{ln{ffj)/dz^)z=o] , equals 
a for a Gaussian, and is infinite for the square function. 
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Our calculations show that, for the case of the square 
influence function (inflnite E), patterns of nontrivial am- 
plitude appear for all values of the cut-off interval of the 
square, the amplitude of the peaks being not uniform but 
presenting a periodic modulation. For a Gaussian influ- 
ence function (flnite E), the patterns exhibit a curious 
feature in the periodic boundary condition case. Two 
critical values of the width E are seen separating trivial 
patterns with vanishing amplitude from nontrivial pat- 
terns with substantial amplitude This is displayed 
in Figs. Q and [21 The critical width depends linearly 
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FIG. 2: Dependence of the amplitude of the patterns on the 
ratio of the width of the influence function to the domain size. 
Inset, which depicts the small S/L behavior on an expanded 
scale, shows an apparent transition around the value 2/9 of 
E/L. The arrows mark the three patterns depicted in Fig. 1. 
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FIG. 1: Dramatic difference in the amplitude of the steady 
state patterns with Gaussian influence function and periodic 
boundary conditions, for E/L = a) 0.205, b)0.245 and c)12. 
Here, u{x) is plotted in units of a/b, and x in units of L. 

on the domain size L (with a slight deviation for small 
domains) making it possible to display our results with 
E/L on the x-axis. The results displayed are for a system 
size L of 100 sites with D = 1 x 10~^, a = 1, 6 = 1, in 
arbitrary but consistent units. A plot of the pattern am- 
plitude Aq (see Fig. ^ shows a striking transition around 
the value 2/9 of the ratio E/L. 

Is this sharp difference between square and Gaussian 
influence functions the result of the cut-off inherent in 
the former? In particular, could the sudden rise of the 
pattern amphtude in the inset of Fig. 2 occur because of 
the onset of the natural cut-off imposed on the Gaussian 
by the flnite domain size? In order to answer these ques- 
tions, we considered domains large enough (so that the 
domain size was unimportant) and used a combination of 
a Gaussian and a square influence function, i.e., a Gaus- 
sian of width E, (equivalently, range cr), multiplied by a 
pulse so that it vanishes abruptly beyond a cut-off Xc- 
We found results which suggest that the non-negligible 
patterns are indeed associated with the cut-off nature of 



a square function or the cut-off imposed on a Gaussian 
by the finite size of the domain: large-amplitude patterns 
appeared for cut-off Gaussian for any value of E provided 
Xc was small enough. 

In order to understand the interplay of E and Xc bet- 
ter, we vary them independently and construct a phase 
plane (E,Xc) in Fig. 01 The separation of the region 
of patterns with large amplitude (displayed shadowed) 
from that of patterns of negligible amplitude is clear and 
occurs along a line of slope 9/2 [l3|. This value of the 
slope is significant in fight of the fact that in Fig. 2, the 
transition is seen at E/L f« 2/9. This clarifies that it is 
the finite domain size that imposes a cut-off on the Gaus- 
sians considered in Fig. 2, and that the domain size L 
there plays the precise role of 2xc in Fig. 3. 

The cut-off Gaussians or the square infiuence functions 
possess an abruptness feature which would not be present 
in a physical system. To ensure smoothness, we borrow 
from [131 and consider an infiuence function (for periodic 
boundary conditions) 

Zr g[r,a) 

_ l ^g{r,a) r(l/r + l) 
" V r r(l/r + 3/2)' 
g{r,a) = (2 + 3r)cr^ (6) 

Here r is non-negative, T is the Gamma function, and 
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FIG. 3: Separation of large-amplitude patterns (shadowed re- 
gion) from negligible-amplitude patterns for the cut-off Gaus- 
sian influence functions with width S and cut-off length Xc- 
The domain size is L = 100. Units for S and Xc are arbitrary 
but identical. 



the influence function h cut-off at with 
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FIG. 4: Counterpart of Fig. 3 for the influence function of 
Eq. Inset shows /cr.r in the periodic boundary condition 
case for several values of r = 100 (dashed), 1 (dotted), and 
.01 (full); respectively, S = 70, 7, and .7. The square and the 
Gaussian emerge as particular cases. Units as in Fig. 3. 

The limit r = is the Gaussian. Note that denotes 
a natural cut-off of this function while Xr in the Gaussian 



case is imposed externally by multiplying by a symmetric 
square function of width 2xc. The width E, the range a, 
and the parameter r are related through 
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as in ref. [12], but analysis with 
's does not add anything important 



r(-i/r) 
such long-tailed f^^ 
to our present study, the results being similar to those 
for Gaussians. In addition to the smoothness property of 
these functions, they have the feature that they reduce 
to a square or a Gaussian in the respective limits r — > 
and r — > oo (see inset in Fig. It is therefore pos- 
sible with their help to study for a given cut-off length 
the effect of varying the width (central curvature) and 
vice-versa. We fix a value of the natural cut-off well 
within the domain size L, and plot in Fig. 4 the counter- 
part of Fig. 3 for the class of influence functions given by 
Eq. ©. The two plots are similar in that the regions of 
large-amplitude and small-amplitude patterns are sepa- 
rated cleanly by what appears to be a straight line in the 
phase space. The slopes are different (9/2 for Fig. 3 but 
about 8 for Fig. 4). We conclude that the cut-off length 
and the central curvature are both relevant to the forma- 
tion of patterns in an interrelated manner shown by the 
phase plots. It should be noticed that while the Gaus- 
sian has an inflnite cut-off Xc and a flnite width S, the 
square has a finite cut-off Xc and an infinite width (cur- 
vature at center). Our introduction of Gaussians with 
an external cut-off provided in the first part of our inves- 
tigation Gaussians wherein both control quantities were 
finite. Our introduction of the influence functions de- 
flned in Eq. © similarly provided square-like functions 
wherein both control quantities were finite. The inset of 
Fig. 4 makes it particularly clear that we can produce 
influence functions for this case with a fixed cut-off and 
S ranging from to oo. 

We have performed a number of further studies of pat- 
tern formation from the long-range Fisher equation (2). 
They include extension to 2-dimensional systems, typi- 
fled by Fig. [HI in which patterns are shown for periodic 
boundary conditions and a cut-off function of the kind 
described in Eq. JHJ. We have also studied the effects of 
boundary conditions. As an example we display in Fig. 
El the lower-symmetry patterns for zero-flux conditions 
that appear to have additional modulation in the peaks. 
We have also investigated the time evolution of the pat- 
terns and found it to have considerable complexity even 
for simple initial conditions. We will report these various 
results elsewhere. Here we would like to draw attention 
to the features of the influence function, cut-off length 
and central curvature or width, whose ratio appears to 
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aspects of the patterns. We hope that future investiga- 
tions will clarify at a basic level the picture presented 
by our numerical findings and make possible predictable 
manipulation of patterns in real systems such as bacterial 
colonies by controlling the influence function via control 
of the flow of nutrients and/or chemotactic substances. 
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FIG. 5: Typical steady-state pattern in a 2D system with 
periodic boundary conditions and a cut-off influence function. 
Lighter areas correspond to larger values of u. Parameters are 
arbitrary. 



determine whether the steady-state patterns have negli- 
gible or sizeable amplitude, the critical value of that ratio 
being different for different families of influence functions. 
We have systematized our numerical findings through the 
phase plots in Figs. ElandEl Our use of simple infiuence 
functions, such as the square, the Gaussian, and the cut- 
off Gaussian, as well as of richer functions as in ref. 
which reduce to these forms, has resulted in our being 
able to focus on what features are responsible for what 




This work is supported in part by the Los Alamos Na- 
tional Laboratory via a grant made to the University 
of New Mexico (Consortium of The Americas for In- 
terdisciplinary Science) and by National Science Foun- 
dation's Division of Materials Research via grant No 
DMR0097204. 



FIG. 6: Typical steady state patterns with cut-off influence 
function and zero-flux boundary conditions. 
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